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ABSTRACT 

We study the environmental dependence of stellar population properties at 
z ~ 1.3. We derive galaxy properties (stellar masses, ages and star formation his- 
tories) for samples of massive, red, passive early-type galaxies in two high-redshift 
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clusters, RXJ0849+4452 and RXJ0848+4453 (with redshifts of z = 1.26 and 1.27, 
respectively), and compare them with those measured for the RDCS1252. 9-2927 
cluster at z=1.24 and with those measured for a similarly mass-selected sample 
of field contemporaries drawn from the GOODS-South Field. Robust estimates 
of the aforementioned parameters have been obtained by comparing a large grid 
of composite stellar population models with extensive 8-10 band photometric 
coverage, from the rest-frame far-ultraviolet to the infrared. We find no varia- 
tions of the overall stellar population properties among the different samples of 
cluster early-type galaxies. However, when comparing cluster versus field stellar 
population properties we find that, even if the (star formation weighted) ages 
are similar and depend only on galaxy mass, the ones in the field do employ 
longer timescales to assemble their final mass. We find that, approximately 1 
Gyr after the onset of star formation, the majority (75%) of cluster galaxies have 
already assembled most (> 80%) of their final mass, while, by the same time, 
fewer (35%) field ETGs have. Thus we conclude that while galaxy mass regulates 
the timing of galaxy formation, the environment regulates the timescale of their 
star formation histories. 

Subject headings: galaxies: clusters: individual: RXJ0849-F4452, RXJ0848-F4453, 
RDCS1252. 9-2927 — galaxies: high-redshift — galaxies: fundamental parame- 
ters — galaxies: evolution — galaxies: formation — galaxies: elliptical — cos- 
mology: observations 



Introduction 



Galaxies reside in environments that span a wide range of density. In order to understand 
the physical processes that drive their evolution is important to test for systematic differences 
between galaxies in various environments. Many authors have shown that population density 
plays an important role in determi ning many galaxy prop e rties, such as star formation 
rate, gas content and morphology (iKodama fc Bowerl l200ll : iBalogh et al.ll2002[ ). Several 
mechanisms have been proposed by t heorists to account fo r these effects, such as ram pressure 
stripping, rnerger s and tidal effects (iGunn fc GottI Il972l : iDressler et al.l 119971 : iMoore et al. 
19961 . ll998[E999l ). 



More than half of all stars in the local Universe are found in massive spheroids (e.g.. 
Bell et al.l ( 120031 )). The history of mass assembly and star formation of these Early- Type 
Galaxies (ETGs) are among the most actively pursued elem ents in galaxy evolution studies, 
and form the basis for models of massive galaxy formation (jRenzinill2006l ). 
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Galaxy clusters provide an ideal laboratory to study ETGs, the dominant galaxy population 
in clusters, even beyond redsliift one. In addition, large field surveys have multiplied in 
recent years, extending the baseline over which environmental effects can be studied . 

Studies indicat e that the most massive ETGs in the f i eld may be amongst the oldest 
at any given epoch ( Cimatti et al.ll2004 : 



2005: Juneau et al. 



Fontana et al. 



20051 : Idi Serego Alighieri et al. 



2004 



2005 



Pannella et al. 



Saracco et al. ||20p4j: iTreu et al. 

2009h . They have 



evolved mainly passively since I and, at least the most massive ones, as slowly as cluster 



2000; van Dokkum &: Franx 



galaxies (Ivan Dokkum fc Franx 



199 



20011 : 



Bernardi et al. 



Treu et al.ll2001 



1998 



Treu et al. 



19991 



van Dokkum fc Stanford 



Kochanek et al. 



20031 : iDe Propris et al. 



20071 ). Similarly, tight constraints have been placed on t he scatter and s 



magnitude relation (CMR), the Fundamental Plane (F 



2005 



van der Wei et al.l 120051 : Idi Serego Alighieri et a. 



20071 ) and Balmer line strengths ( IClemens et al. 



2006 



) (|j0rgensen et al. 



2006! 



ope of the color- 



2006 



Treu et al. 



van Dokkum &: van der Marel 



Sanchez-Blazquez et al. 20061 ) of the 



ETGs populati ng massive clusters, i ndicating a very high formation redshift for the ETGs in 



clusters as well ( 


Rettura et al. 


2006 


Eisen 


lardt et al.l 


2008; 


Rettura et al. 


2010; 


Gobat et al. 


2008; 


Mei et al. 


2009; 


Collins et al. 


2009 




Rosati et al 


2009 


)• 



The most distant clusters 



dictions (jPapovich et al.ll2010 



ETGs at 0.8 < z < 1.5 fBlakeslee et al 



Lidman et al. 2008; Hilton et al. 



mown to date provi de the strongest leverage on model pre- 



Tanaka et a. 



2009 



2010af). Specificallv . ther e is a tight CMR of 



20031: iBlakeslee et al. 2006; Mei et al. 2006a b 



StrazzuUo et al.ll2010[ ). a slowly evo lving -R'-band lu- 



2004 



minosity function see mingly at odds with hierarchical merging scenarios (iToft et al 
Strazzullo et al. 2006), and a tight and slowly evolving FP out to ^ ~ 1.2 (jHolden et al. 
20051 ). The existence of such massive, passively evolving galaxies already at z > 0.8 



needs to be reconciled with the evidence for massive (often dusty) st ar-forming popula- 



tions found at z ^. 2.2 both in the field and in overdense environments (ISteidel et al 



Adelberger et 


al. 


2005 


Tanaka et al. 


2010bh. 



2005 



2009 



Previous studies based on optical and infrared data have placed some constraints on the 
star formation histories of ETGs in distant clusters based on the scatter in their CMRs: the 



best-fit models yield formation epochs of 2 < 2;^- < 7 (IBlakeslee et al.ll2003l ). but a greater 
accuracy could not be achieved based on the available rest-frame optical and infrared data 
alone. The reason for this is that age-dependent features in galaxy spectra (e.g., the 4000 A 
break or th e 1.6 fj,m bump) are too broad or evolve too slowly to be used as adequate ag e 



indicators (IBurstein et al.l 1 19881 : iTantalo et al.l Il996l : iMarastonI l2005l : iRettura et al.l I2OIOI ) . 

In order to better constrain the formation epoch [zj) of ETGs we need to probe the much 
sharper and fast-evolving rest frame UV-optical colors that will allow us to determine Zf with 
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much greater accuracy than hitherto achie v ed, en abling us to better constrain the entire SFH 
of ETg£| (see also Fig. 2 of iRettura et~aD toid )) 



In order to provid e a key test of the paradigm of an accelera ted evolution in the highest 
density environments ( jPiaferio et al.ll200ll : iThomas et al.ll2005l ). in this work, on the basis 
of photometric data available over the rest-frame wavelength range 0.15 — 2yum, we compare 
stellar masses, ages and inferred star formation histories of ETGs found in three massive X- 
ray detected clusters at ^ ~ 1.3 (RXJ0849+4452, RXJ0848+4453, RDCS1252.9-2927) with 
a sample of field contemporaries drawn from the GOODS South Field. 

The structure of this paper is as follows. The description of our datasets, cataloging 
and sample selection is described in §2. In §3 we describe our methods in inferring star 
formation histories, ages and masses from stellar population analysis. The results of our 
study are discussed in §4, while in § 5 we summarize our conclusions. 



We assume a flf 



0.73, Qr 



0.27 and Hq = 71 km • s ^ • Mpc ^ flat universe 



( ISpergel et al.l 120031 ) . and use magnitudes in the AB system throughout this work. 



2. Description of the data 

This work is based on data collection programs on the three galaxy clusters at z ~ 1.3 
which are amongst the most extensive spectroscopic and photometric surveys available over 
the wavelength range 0.35 — 4.5;um. To compare cluster galaxy properties with the field, we 
also take advantage of similar- quality archival data available for the GOODS-South Field. 
Tab. [T] summarizes the rich and homogeneous dataset employed in this study. 



2.1. RX J0849+4452 and RX J0848+4453 



The Lynx Supercluster, is the highest redshift supercluster known today (INakata et al. 



20051 ) , with two central X-ray detected galaxy clusters and three spectroscopically confirmed 
surrounding groups. This work focuses in particular on the cores of the two main clusters, RX 
J0849+4452 (hereafter. Lynx E) a nd RX J0848+4453 (hereafter. Lynx W), both detected in 
the ROS AT Deep Cluster Survey (IRosati et al.lll999[) and spec troscopicallv confirmed at z 
= 1.261 (IRosati et al.lll999[ ) and z = 1.273 ([Stanford et al.lll997l ). respectively. Furthermore, 



^Note that by z >, 1.3 not enough cosmic time has elapsed for HB stars to produce the "UV upturn' 
|Yi et aLlll99"7[ll99<£ 
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in two accompanying papers (IRaichoor et al.ll201ll : iMei et al.ll201l[ ) we also compare ETGs 
stellar population properties in the two cluster cores with those of the ones in the surrounding 
groups. 

Both Lynx E and Lynx W were observed in the u' filter with the blue channel of 
LRIS on Keck I in two separate epochs. On March 4th, 2003, exposures totaling 260 min 
were obtained in conditions that were photometric with 0.8" seeing. A further 240 min was 
obtained during the night of March 1st, 2008, with clear skies and a typical seeing of 0.9" . The 
two epochs of data were reduced separately using standard IRAF tasks, where we subtracted 
the bias level, flat-fielded each exposure using a series of twilight flats obtained on each night, 
astrometrically calibrated and aligned each image by comparison to star positions from the 
HST GSC2.2 catalog and then refined with a larger number of stellar positions from the 
SDSS DR6 catalog. Finally, the images were fiux-calibrated, also correcting for atmospheric 
extinction, using exposures of various photometric standard stars obtained across a range 
of airmasses each night, and then combined. The two sets of images were taken at different 
position angles, and so the combined frame has an effective exposure time of 8 hrs 20 min 
across much of the field, but a shallower depth at the edges of the field where the frames do 
not overlap. The effective resolution in the combined ultraviolet image is 0.95". 

In order to sample the optical-to-infrared wavelength domain, we have obtained data 
in the following passbands, R, imbw^ zfs5olp, J, Ks, and Spitzer /IRAC channels, 3.6/im 
(chl) and 4.5/im (ch2). 



We refer the reader to the companion paper by (IRaichoor et al.ll201lh for a more detailed 
description of the optical and infrared dataset used in this work (see also their Tab.l). 

Here we remind that the R-band images were obtained with the Keck/LRIS instrument. 

The 

^F775H/ and 2;i;'850Lp-band were observed in the F775W and F850LP filters provided by 
the ACS- Wide Field Camera (WFC) on board of the Hubble Space Telescope. We note 
that these two observing filters have been purposely chosen to bracket the 4000A-break of a 
model elliptical galaxy at z = 1.26 — 1.27 f lMei et al.ll2006al Jbl). 

The near-infrared imaging data (J, i^'s-band) were acquired with the FLAMINGOS instru- 
ment available at the Kitt Peak National Observatory (KPNO). 

Mid-Infrared imaging was obtained (in two channels at [3.6yum] and [4.5^m]) with the IRAC 
camera on board of the Spitzer space telescope. 
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2.2. RDCS1252. 9-2927 and the GOODS-South field. 



This work builds on datasets and analyses already performed on ETGs belonging to 
the X-ray luminous cluster, RDCS1252.9-2927 (hereafter, CL1252), and to the GOODS- 
South field, centered on th e so-ca l led Chandra Deep F ield South field (hereafter, CDFS), as 
reported bv lRettura et al.l fcoioh : ICobat et all tooi ). 



We refer to the aforementioned papers for more details in the data reduction. Here 
we note that the data we have employed for the CL1252 cluster consist of deep imag- 
i ng in 10 bands : FL T/VIMOS (f/), VLT /F 0RS2 (B, V, ffl, E ST /ACS (iF775iy,^F850/p) 



(iBlakeslee et al.l 120031 ). VLT /ISAAC {Js,K,) flLidman et al.ll2004f ). Spitzer/I KAC (3.6Atm, 
4.5/xm ), as well as spectroscopic data taken with VXr/F0RS2 and published in lDemarco et al, 
J2007h . 



The archival data for the comparison fi eld, the CDFS, comprises deep imaging in 9 
bands: FLr/VIMOS j U) jNonino et al. l l2009l ), i/5r/ACS {B^^^,,w , VpmQW, iri^w, zfs^olp) 
Jciavalisco et al.lbooi . FLT/ISAAC, {J,K,) jRetzlaff et al.lboioh . Spitzer/IRAC (3.6/im, 
4.5/zm)[§, as well as spectroscopic data taken with the VLr/F0RS2 30 01 grism by the ESO- 
GOODS surve\R JVanzeUa et al.l boosl booei boosi : ISalestra et aDboioh and the K20 survey 
(Cimatti et al. 2002). 



2.3. Cataloging of observations and samples selection 



The resulting datasets for the three clusters and the field have homogeneous depths and 
wavelength coverage, allowing the application of similar selection criteria for both samples. 
This is a crucial point for stellar population studies, so that we do not have different levels 
of systematic biases in the analyses. 

The data allow the reconstruction of galaxy spectral energy distributions (SEDs) by 
entirely sampling the relevant spectrum range emitted by all the different stellar populations. 
For each observing band we produce catalogs from matched aperture photometry, with 
an aperture o f radius 1.5" and an apertur e correction out to 7" radius, as described in 
Rettura et all fcood ): iRaichoor et al.l 3201l[ ). 



^CDFS imaging is publicly available through the GOODS collaboration web-site: 
http: / / www.stsci.edu/science/goods/ , 

^Spectroscopic data are publicly available through the web-site: 

http://www.eso. org/sci/activities/projects/goods/MasterSpectroscopy.htm^. 
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The availability of 8 to 10 passba nds spanning such a. large wavelength range enables 
the estimate o f accurate stellar ma sses (iRettura et al.ll2006l ). ages and star formation histo- 
ries of ETGs (IRettura et al.l 120101 ) and enables us to directly compare galaxy properties of 
homogeneously selected samples of ETGs in both environments. 

Specifically, we measure stellar population parameters of Lynx E and Lynx W early- type 
cluster galaxies and compare them with those found for similarly-selected samples of cluster 
and field contemporaries drawn from CL1252 at z = 1.237 and CDFS aX z = 1.237 ± 0.15. 

We estimated stellar masses from both Ks band observations and SED fitting . The two 
estimates give similar mass limits. The depth of the FLT/ISAAC and ifPA^O/FLAMINGOS 
i^'s-band images and the extended multi- wavelength data for all fiel ds allows us, i n fact , 
to define comp l ete m ass-selected samples. The reader is referred to iGobat et al.l ( 120081 ): 



Raichoor et al.l (120111 ) for more details. Here, we note that overall photometric completeness 
is obtained if we limit our analysis to stellar masses larger than Mum = 5 ■ 10^°M(rO. 

A selection of Lynx E, Lynx W and CL1252 passive ETGs along th e cluster red sequence 



Mei et al. 


2006a 


b, 


2009. 


2011) 



zf&^olp > 0.8 (IBlakeslee et al.l 12003 



For CL1252, in the spectroscopic sample of iDemarco et al.l (120071 ). there are 22 red-sequence 
galaxies {ipn^w — zfs50lp > 0.8) with > Mum, of which 18 are classified as passive 



Lynx E and Lynx W, in the combined photometric and spectroscopic samples ( 


Mei et al. 


2006a 


:lstanford et al. 


1997; 


Rosati et al. 


1999; 


Stanford et al. 


2001 


van Dokkum & Stanford 



2003 



Holden et al.l 120111 ) there are 21 red-sequence galaxies with M* > Mum, of which 17 
are classified as passive ETGs. Specifically, there are 10 passive spectroscopically-confirmed 
ETGs in Lynx E, and 3 in Lynx W. We remark that, to improve the statistics, we have 
also included into the Lynx ETGs sample, 3 photometrically-confirmed members of Lynx E 
and 1 photometrically-confirmed member of Lynx W both with iprrsw — zpsboip > 0.8 and 
M^, > Miim, thus resulting in a total sample of 17 Lynx E+W ETGs. It is important to 
note that the removal from our sample of these 4 galaxies, whose cluster-membership was 
assessed on the basis of a (8-9 band) photometric-redshift only, would not change the main 
conclusions of this study. 

For the corresponding CDFS sample of field contemporaries, the same criteria yield 27 passive 
ETGs in CDFS with F0RS2 spectra giving redshift in the range z = 1.237 ± 0.15 (for more 



Assuming Salpeter Initial Mass Functio n (IMF), with l ower a nd upper cutoffs , mL = O.IMq and mu = 
100 Af©, respectively. We refer the reader to IRettura et al.l (|2006l ): iRaichoor et al.l (j201lt ) for a discussion on 
the dependence of inferred stellar population properties on the actual choice of the IMF. 
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details see also ICobat et al.l (120081 ): iRettura et al.l (120101 )1 

We also remark that, adopting the classification scheme of Postman et al. (2005), visual 
morphological analysis of the HST/ACS images available for all samples indicates that all 
our red, passive ETGs also show typical elliptical or lenticular morphology. 

As reported in Rettura et al. 2010, we also remind that the spectroscopic follow-up for 
CL1252 is more complete at the low-mass end than in CDFS. Thus, our sample of ETGs in 
CDFS is likely to be more incomplete at the low-mass end than the CL1252 and Lynx ones. 
We will return to this point when discussing our results in Section 4. 



3. Data Analysis: derivation of stellar population properties 

The data described in the previous section are used to infer fundamental physical prop- 
erties of similarly mass-selected samples of early-type galaxies (13 ETGs in th e Lynx E, 4 in 



Lynx W, 18 in CL1252 and 27 in the CDFS) . Adopting a similar approach to iRettura et al. 



(I2OO6I . I2OIOI ). we derive stellar masses, star formation histories and ages for each ETG in Lynx 
E and Lynx W using multi- wavelength PSF-matched aperture photometry from 8 passbands, 
from observed w'-band to observed A.bfim. For each galaxy, we compare the observed SED 
with a set of comp o site s tellar populations (hereafter, CSP) t emplates compu ted with the 
Bruzual fc CharlotI (120031 ) models, assuming solar metallicity, ISalpeteii (119551 ) Initial Mass 
Function (hereafter, IMF) and no dust. 



Similarly to IRettura et al.l ( l2006l ) we checked the effect of dust extinction on the best-fit 
stellar masses by inc luding a fourth free parameter, 0.0 < E{B — V^) < 0.4, following the 
Cardelli et al.l (Il989[ ) prescription. Performing the SED fit on 13 Lynx E+W ETGs (all 
with complete u'-to-i.dfim coverage) and on 22 CDFS ETGs (all with complete U-to-i.bfim 
coverage), we find that in ~ 50% of the cases E{B — V) < 0.05 gives the best fit. In the 
remaining cases values of E{B — V) < 0.2 are found and no particular trend with mass nor 
environment is found either, hence supporting the validity of the dust-free assumption we 
make throughout this work. 

For our CSP models, we assume the following grid of exponentially-declining star for- 
mation history (SFH) scenarios, \E'(t, r): 



Mr. 







where 0.05 < r < 5 Gyr, SFRq is the initial SFR at the onset of star formation, and 



- 9 - 



t is the time since the onset of star formatioiifl of the stellar population model formed at a 
lookback time T(z) + 1, with, e.g., T 8.67 Gyr at the epoch of observation [z ^ 1.3). 

In determining galaxy model ages, masses and star formation histories from SED fitting, 
is important to understand how much our estimates could possibly be af f ected by "age- 
metallicity"!^ and "age-SFH" degeneracies. We note that in Rettura et al. ( 2010 ) we have 
demonstrated that the use of photometric information coming from the rest-frame UV is 
crucial to distinguish the different parameters of the stellar population modeling, allowing 
us to break the "age-SFH degeneracy" . It is also important to note that the rest-frame UV 
remains also largely un affected by the "age-metallicity" degeneracy, which plagues optical 
studies ( )Worthevlll994j ). 



Furthermore, the reliability, at z > 1, of stellar popula tion parameters inferred in the 
rest-frame NIR regime {Xobs ~ 2/im) has long been debated (lMarastonlll998[ ) because of the 
different implementation of the relevant short- duration thermally p ulsating (TP) AGB p hase 
in the different stellar population synthesis models. However, in iRettura et al.l (120061 ) we 
have shown that different stellar population models actually yield overall consistent stellar 
masses, within the typical errors, for ETGs at 2; ~ 1.0. It is clear that the inferred SFH 
(and age) is model-dependent and is generally not a unique solution, howev er the relative 



differe nce in the underlying stellar population is still significant. In fact, in iRettura et al. 



( 120101 ) we have shown that when comparing cluster vs. field rs and ages, the relative dif- 
ferences among samples (or lack thereof, as in their Fig. 8) were similar regardless of the 
actual stellar population code used, not affecting their conclusions. 

Hence, since our study is indeed aimed at constraining the relative difference of galaxy stellar 
population parameters in a Cluster vs. Cluster and Cluste r vs. F ield fashion, we will only 
show here the results obtained with the iBruzual fc CharlotI (120031 ) models, i n order to facil 



itate direct comparison with previous studies in the literature. Note that in iRaichoor et al. 



( I2OIII ) we also present an extensive discussion on systematics due to choice of the stellar 
population models at z ~ 1.3 . 

To account for the average age of the bulk of the stars in a galaxy, we refer throughout 
this paper to star-formation weighted ages, < t >sfh, defined as: 



< t >SFH^ 



J^{t-t'Mt\r)dt' 

j'^it',T)dt' 



(2) 



^The range of acceptable is for a given galaxy has been limited by the age of the universe at its observed 
redshift. 

^We employ the working assumption that the most-massive ETGs have all solar metallicities. 
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Assuming \l/(t',r) as in Eq. [T]we obtain, 



< t >SFH- 



T + T ■ e 



(3) 



where t is the time elapsed since the onset of star formation (SF). 

By comparing each observed SED with these atlases of synthetic spectra, we construct 
a 3D space spanning a wide range of star formation histories, model ages and stellar 
masses. The galaxy mass in stars M*, the inferred < t >sfh and the r of the models giving 
the lowest are taken as the best-fit estimates of the galaxy stellar mass, age, and SFH 
timescale throughout this work. We note that this procedure resu lts in typical errors for 
galaxy ages of ~ 0.5 Gyr, and for r of ~ 0.2 Gyr ( iRettura et al.ll2010l). Tvpical un certainties 
on the mass determination are about ~ 40% (i.e., 0.15 dex) ( iRettura et al.l 120061 1. 



4. Results and Discussion 

In the top panel of Fig. [1] we plot the inferred stellar masses as a function of the 
U — B rest-frame color for the mass-selected samples of CL1252, Lynx E and Lynx W 
ETGs, resulting in a very similar distribution. In the bottom panel of Fig. [1] we show the 
rest-frame U — B color - mass diagram of the combined samples of cluster (red circles) 
and field (blue circles) ETGs. To have a very simple estimate of the scatter of the color- 
mass relation, as de rived from the SE P fitting, we derive the scatter around the fit using 



a Tukey's biweigth ( jPress et al.lll992l ). We simply calculated the uncertainty in the scatter 
estimation by bootstrapping on 1,000 simulations. For the Lynx cluster sample, CL1252, 
and field galaxies, we obtain a scatter of 0.031 ± 0.010, 0.042 ± 0.010 ,and 0.050 ± 0.008, 
respectively. These are intrinsic scatters as predicted from the best SED fitting, and do not 
take into account uncertainties in the fitting methods. Field ETGs galaxies are distributed 
around the cluster red-sequence, although they seem to show a larger scatter, especially at 
lower masses. We will need more statistics and better estimates of uncertainties to draw 
a stringent conclusion, since the cluster and field overall scatters could still be consistent 
within the uncertainties. 

As we apply the method described in section §3, we are able to directly compare the 
relative distribution of galaxy stellar population properties for the different samples. 
Stellar masses ranges from to 4 • IO^^Mq except for three CL1252 bright galaxies and 
one galaxy in the field (bottom panel of Fig. [T]). We do not find any of the Lynx clusters 
ETGs to be more massive than 3 ■ IO^^Mq (top panel of Fig. [1]). 



We want to tests the hypothesis that two samples stellar-mass distributions have the 
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same median against the hypothesis that they differ at the 5% significance level. The Mann- 
Whitney U-Test suggests that the distributions are likely to be similar (median probability 
of 44%). We also study the Kolmogorov-Smirnov statistic and associated probability that 
the two arrays of data are drawn from the same distribution, yielding a similar result. A 
high value (67%) of probability shows that the cumulative distribution function of the cluster 
stellar masses is not significantly different from the field one. 

As shown in Fig. [21 we also find that all clusters show the same relative distributions 
of lookback time since the onset of star formation (T{z) + t, top panel) and stellar mass 
assembly timescales, (r; bottom panel). In fact both the Mann- Whitney U-Test (median 
probability of 10%, 21 %) and the K-S statistics (75%, 91% of probability) also elucidate 
that the distributions are very similar. 

This lack of cluster-to-cluster variation may well indicate that these three clusters ex- 
perienced similar evolutionary processes, resulting in similar stellar population properties 
shown by their member ETGs. Since these samples of cluster ETGs are found to be very 
similar from a stellar population point-of-view, we can combine them into one sample of 35 
ETGs, representative of the Cluster environment at z ~ 1.3, apt for comparison with the 
sample of 27 CDFS ETGs, representative of the Field environment. 

Thus we compare the distribution of the star-formation weighted ages (< t >sfh, top 
panel of Fig. [3]) as a function of stellar mass in both environments. We find the distributions 
to be overall very similar, as also shown by the histograms of lookback times since the onset of 
star formation (top panel of FigH]) and by the histograms of the lookback times to < t >sfh 
(bottom panel of FigH]). Note that Mann- Whitney U-Test returns a median probability of 
18% and the K-S statistics a 60% of probability that the last two data values are drawn 
from the same distribution. These results imply that no significant delay in relative age is 
found for ETGs in either environments. Using BC03 models, we find that ~ 80% of ETGs 
have SF-weighted ages in the range 3. 5 =fc 1.0 Gyr in both cluster and field, in qualitative 
agreement with previous CMR studies (IBlakeslee et al.ll2003t iMei et al.ll2006al Jbl). 



In the bottom panel of Fig. |3] we show, for both samples, the dependence of star- 
formation weighted ages, < t >sfh, on stellar mass. It is evident that the ages of ETGs 
only depends on their mass, i.e. the halo mass in which t hey reside, which i s in agreement 
with the so-called downsizing scenario of galaxy formation (jCowie et al.lll996l ). We note that 
this scenario have also been reconciled, with the introduction of various forms of so-called 
feedback m echanisms, with the r ecent versions of s emi-analytic mode ls based on AC DM 
cosmogony JPe Lucia et~aDl2006l : ISower et al.lbood : iMenci et al.lboosh . 



To summarize, we find that cluster galaxy (star formation weighted) ages have the same 
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relative distribution as their field contemporaries: no significant delay in their formation 
epochs is found within the errors (~ 0.5 Gyr). We remind that this result is in rem ark- 



ably good agreement wi t h the ones derived by Ivan Dokkum fc van der Marell (120071 ) and 



di Serego Alighieri et al.l (120061 ) from the evolution of the M/ L ratio. However our method 
is able to extend this kind of analysis to less massive system than the typical targets of FP 
studies at z ~ 1. 

Despite of the fact that cluster and field galaxy formation epochs are found to be 
similar, it could still be possible that the timescales of their SFH are significantly different. 
As mentioned above, the data have shown that the distributions of cluster and field optical 
colors were slightly different. As a function of the stellar mass, cluster galaxies are found to 
lie on a very tight red-sequence, while those in the field populate the same color-sequence 
with a slightly larger scatter (Fig. [I]). 

This piece of evidence finds a natural explanation in the framework of our modeling. As 
shown in the top-left panel of Fig. [5l as a function of stellar mass, we find that field ETGs 
show longer SFH timescales than their cluster contemporaries, which tend to assemble their 
mass with the shortest r at any given mass. In fact the Mann- Whitney U-Test confirms 
that the distributions are very likely to be different (median probability of 0.07%). The 
Kolmogorov-Smirnov statistic also yields an associated probability that the two distribution 
of rs are drawn from the same distribution of just 0.8%. 

If galaxies continue to form stars following the exponential decay of Eq. (H we can derive 
the time, to.s needed to form 80% of the stars they would have at 2; = 0. As shown by top- 
right panel of Fig. [5] and bottom panel of Fig. [6]), we find that ~ 1 Gyr after the onset of 
star formation 75% of cluster ETGs have already assembled most of their of their final mass. 
By the same time, less than 35% of field ETGs have. 

The bottom-left panel of Fig. [5] shows the relative distributions of inferred initial star 
formation rate, SFRq, as a function of stellar mass. More than 30% of cluster galaxies 
have higher SFRq and smaller rs than any other field contemporaries, indicating they have 
experienced a much more intense star formation at early times than field galaxies. As a result, 
at the time of observation, cluster ETG specific star-formation rates, defined as the SFR 
divided by the galaxy stellar mass, SSFR, are lower than those of the field ones (bottom- 
right panel of Fig. [5]). Indeed the Mann- Whitney U-Test confirms that the distributions of 
SSFR are likely to be different (median probability of 1.3%). The Kolmogorov-Smirnov 
statistic yields an associated probability of 0.5% that the two distribution of rs are drawn 
from the same parent distribution. 

As summarized by the histograms in Fig. [6l the distributions of r and to.s in clusters 
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and in the field are found to be very different. Tliese results clearly indicate a dependence 
of ETGs star formation histories on the environme nt, in agreenient with previous studies 



based on a compa rison of a single cluster to the field ( IGobat et al.ll2008l : iRettura et al.ll2010 



Menci et al.l 120081 ) . 



As discussed in iRettura et al.l (I2OIOI ) , we recall that our field sample is more deficient 
in lower mass objects than the cluster sample because of a less extensive follow-up of spec- 
troscopic targets in the field as compared to that in the cluster. 

However, even if the field sample were corrected for completeness, this would likely result 
in a larger fraction of field ETGs at low mass, which are the ones that we find with longer 
rs and to.ss, higher SSFRs and lower SFRqs. Hence a correction would actually amplify 
the difference between the typical mass assembly timescales of the two samples, and so not 
affect our conclusions. 



5. Conclusions 

We have studied the environmental dependence of early-type galaxy stellar population 
properties at z ^ 1.3 based on high-quality, multi-wavelength photometry available in 8-10 
passbands from U to [4.5/im]: sampling the entire relevant domain of emission of the different 
stellar populations, from rest-frame far-ultraviolet to the infrared. 

We have similarly analyzed mass-selected samples of ETGs belonging to the cores of 
three massive clusters (RXJ0849+4452 at z=1.261, RXJ0848+4453 at z=1.273, RDCS1252.9- 
2927 at z=1.237) and compared their stellar population properties with those measured in 
a similarly selected sample of field contemporaries drawn from the GOODS-S survey. 

We have derived stellar masses, (star formation weighted) ages, and star formation 
histories, parameterized as timescales, r, from exponentially declining CSP model templates 
built with BC03 models for samples of massive (M > 5 ■ 10^°Mq), passive (no-emission line 
in their spectra) ETGs at z:i^ 1.3. Apart from a lower level of spectroscopic completeness for 
the least massive field galaxies, that we find not to affect our conclusions, our sample has 
the advantage of being photometrically complete at our mass limit and having galaxy types 
(passive or star-forming members) assigned spectroscopically. 



This wo r k ext ends the analyses performed on a single-cluster by IRettura et al.l (I2OIOI ): 



Gobat et al.l (120081 ) to two more clusters at z ^ 1.3: doubling the size of the cluster ETG 
sample. 

Hence, on the basis of the new data we are able to affirm that: 
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• we find no significant difference in the derived stellar population properties of z ~ 1.3 
ETGs belonging to different X-ray luminous clusters (RXJ0849+4452 & RXJ0848+4453, 
RDCS1252.9-2927). 



Moreover the comparative statistical analysis performed in this work on the ext ended 
dataset corroborates the results obtained in our previous study (IRettura et al.l 120101 ). that 
we summarize as follows: 



• Field ETGs galaxies are distributed around the cluster ETG red-sequence, although 
they seem to show a larger scatter, especially at lower masses. We will need more 
statistics and better estimates of uncertainties to draw a more stringent conclusions, 
since the cluster and field overall scatters could still be consistent within the uncer- 
tainties. However, we remind that the small scatter in cluster environments remains a 



challenge for semi-analytical galaxy evolution models (IMenci et al.ll2008l ). 



• We find no significant delay in the star formation weighted ages of massive ETGs ob- 
served in all clusters and in the field at z~ 1.3. 

The age of ETGs increase with galaxy mass in all environments, which is in agree- 
ment with the downsizing scenario. The site of active star formation must have 
shifted from the most massive to the less massive galaxies as a function of the cos- 
mic time. The formation epochs of ETGs only depends on their mass and not on the 
environment they live in. This result is in rema rkably good agreement with those ob - 
tained from the evolution of the M/L ratio (e.g. Ivan Dokkum fc van der Marell (120071 ): 



di Serego Alighieri et al.l ( 120061 )) 



• However, the data show that cluster and field SFHs are significantly different. Field 
ETGs best-fit models span a different range of timescales than their cluster contem- 
poraries, which are formed with the shortest r at any given mass. This result is quan- 
titatively consistent with the predictions of c urrent galaxy form ation models based on 
the latest rendition of semi- analytic models (IMenci et al.ll2008l ). 



• We find that 1 Gyrs after the onset of star formation 75% of cluster ETGs have already 
assembled 80% (or more) of their final mass, while, by the same time, less than 35% 
of field ETGs have. 

• Accordingly, cluster ETGs specific star-formation rates at the time of observation are 
found to be smaller than those of the field ones, implying that the last episode of star 
formation must have happened more recently in the field than in the cluster. 
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While cluster and field galaxy observed at z~ 1.3 form at a similar epoch in a statistical 
sense, a high density environment appear to be able to trigger a much more rapid and 
homogenous mass assembly event for the ETGs, limiting the range of possible star-formation 
processes. In fact, more than 30% of cluster galaxies are found with higher initial SFRq and 
smaller rs than any other field contemporaries, indicating they have experienced much more 
intense star formation at early times than field galaxies. In low density environments, this 
effect must rapidly fade as ETGs display a much broader range of possible star formation 
histories. 
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Table 1: Cluster and Field photometric datasets. 
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Fig. 1. — Top panel: Rest-frame U — B color- mass diagram of mass-selected samples of 
CL1252 (green circles), Lynx E (red circles) and Lynx W (orange circles) passive cluster 
ETGs. Uncertainties in the stellar mass are ~ 0.15 dex. Bottom panel: U — B color - mass 
diagram of the combined samples of Cluster (red circles) and Field (blue circles) ETGs. Field 
ETGs galaxies are distributed around the cluster red-sequence, although are found with a 
larger scatter. 
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Fig. 2. — Top Panel: Histograms of the lookback times to the onsets of star formation for the 
cluster samples. Uncertainties in ages are ~ 0.5 Gyr. No cluster-to-cluster variation of the 
galaxy ages is found. Bottom Panel: Formation timescales, r, histograms for the different 
cluster samples. No cluster-to-cluster variation of the galaxy formation timescales is found. 
Uncertainties in stellar masses and r are ~ 0.15 dex and ~ 0.2 Gyr, respectively. 
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Fig. 3. — Top Panel: The dependence of the star formation weighted ages on the galaxy 
mass and environment. No dependence of ages on the environment is found. Bottom Panel: 
The dependence of the lookback time to star formation weighted ages on the galaxy mass, 
and environment. Uncertainties in stellar masses and ages are ~ 0.15 dex and ~ 0.5 Gyr. 
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Fig. 4. — Histograms of the field (solid blue line) and cluster (solid red line) lookback 
times since the onset of star formation (top panel) and of the star formation weighted ages 
(bottom panel). No dependence of the ETG ages with the environment. Uncertainties in 
age determination are ~ 0.5 Gyr. 
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Fig. 5. — Top-left Panel: Formation timescales, r, of ETG as a function of their stellar 
mass and environment. Top-right Panel: As a function of stellar mass and environment, the 
diagram of the time, to. 85 needed for a given galaxy to form 80% of the stars it will have at 
z = 0. Bottom-left Panel: Dependence of the initial star formation rate, SFRq, of ETGs 
as a function of their stellar mass and environment. Bottom-right Panel: Specific SFRs, 
SSFR, of ETGs as a function of their stellar mass and environment. The mean error in 
stellar age is 0.5 Gyr. Uncertainties in stellar masses, r and to.s, are ~ 0.15 dex, ~ 0.2 and 
~ 0.4 Gyr, respectively. In all panels, as a function of stellar mass, the median values of the 
different cluster (red dashed line) and field (blue dashed line) stellar population parameter 
distributions are also indicated. 
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Fig. 6. — Top Panel: Histogram of the formation timescales, r, of ETGs in both environ- 
ments. Bottom Panel: Histogram of the times, to.S; needed for a given galaxy to form 80% of 
the stars it would have at z = 0, for both environments. Uncertainties in the determination 
of r and to.s are ~ 0.2 and ~ 0.4 Gyr, respectively. 



